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Abstract Three empirical lines of evidence, (Pqcd, pQCD, dA), from RHIC have 
converged and point to the discovery of a strongly coupled Quark Gluon Plasma. 
The evidence includes (1) bulk collective elliptic flow and (2) jet quenching and 
mono-jet production, observed in Au+Au collisions at 200 AGeV, and (3) a crit- 
ical control experiment using D-l-Au at 200 AGeV. 

1. The Theoretical QGP 

The Standard Model of strong interactions predicts the existence of a new 
phase of matter, called a Quark Gluon Plasma (QGP), in which the quark and 
gluon degrees of freedom normally confined within hadrons are mostly liber- 
ated. Lattice QCD calculations show that there is a rapid rise of the entropy 
density, o-(T), of matter when the temperature reaches T Tc ~ 160 MeV. 
Beyond Tc the effective number of degrees of freedom, n(T), saturates near 
the number of quark and gluon helicity states uqcd = 8c x 2^ + | x 3c x Nf x 
2s X 2qq « 37. The entropy density a{T) = dP/dT = (e + P)/T oc n{T)T^ 
approaches the Stefan Boltzmann limit APs b (T) / 3T. The transition region 
is a smooth crossover when dynamical quarks are taken into account, but the 
width of the transition region remains relatively narrow, ATc/Tc ~ 0.1 [1]- 
[6]. 

The rapid rise of the entropy was predicted long before QCD by Hagedorn 
due to the observed exponential rise of the number of hadron resonances [5]. 
However, the saturation of the number of degrees of freedom near uqc d is a 
unique feature of QCD. Even though the entropy density approaches the ideal, 
weakly interacting plasma limit, lattice calculations of correlators show that the 
QGP is far from ideal below 3Tc. The nonideal nature of this strongly coupled 
QGP is also seen from the deviation of the pressure, P{T), and energy density 
e(T) from the Stefan Boltzmann limit as shown in Fig.(l) from [1]. 
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The equation of state of the QGP, Pqcd{T), is the bulk thermodynamic 
property that can be investigated experimentally via "barometric" observables. 
A measure of its stiffness is given by the speed of sound squared, = dP/ de = 
dlogT/dloga = (3 + (ilogn/dlogT)^^ shown in Fig. (2) . Note that drops 
rapidly below 1/3 as the effective number of degrees of freedom drops when T 
approaches Tc. This softening of the QGP equation of state near Tc is a key fea- 
ture can be looked for in the collective hydrodynamic flow patterns produced 
when the plasma expands. 



Figure 1. A recent Lattice QCD calculation [1] of the pressure, P{T) /T^, and a measure of 
the deviation from the ideal Stefan-Boltzmann limit (e(T) - 3P{T))/T^. Note that the scale 
on both graphs has not been corrected for finite lattice volume effects: see [2] for discussion. 

Another distinctive feature of the QGP phase diagram is shown in the right 
panel in Fig.(2). Recent lattice QCD calculations [1] have begun to converge 
on numerical evidence that the QGP may have a second tricritical point [4, 6] 
at moderate baryon densities with fis = 360 MeV and T Tc. 

QCD predictions of the QGP phase date back thirty years [7] and followed 
immediately after the discovery of asymptotic freedom of QCD. The experi- 
mental strategies to search for new forms of dense matter also date back thirty 
years when T.D. Lee proposed "vacuum engineering" [8]. It was then also real- 
ized by W. Greiner and collaborators [9-11] that extended regions of dense nu- 
clear matter can be formed in high energy interactions of heavy nuclei, and that 
the measurement of collective flow patterns will provide the novel barometric 
probes of the equation of state of ultra-dense matter. The hunt for the QGP and 
other phases of nuclear matter has been underway since that time using sev- 
eral generations of higher energy accelerators, BEVALAC, AGS, SPS, and now 
RHIC, and covering an impressive energy range y^— 2m7v = 0.2—200 AGeV. 
In three years, LHC is expected to start vacuum engineering at 5500 AGeV. 
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Figure 2. Important features of the QGP equation of state. The speed of sound [3] = 
de/dP drops below 1/3 for T < 2Tc ~ 300 MeV. Right panel shows a current estimate of the 
location of the tricritical point at finite baryon density [1] 



The first conclusive evidence for (highly dissipative) collective nuclear flow 
was seen at the BEVALAC in 1984 [12], and at AGS and SPS there after. 
However, the first conclusive evidence for nearly dissipation free collective 
flow obeying Pqcd had to await RHIC. In these lectures, the discovery of the 
novel low dissipation elliptic flow pattern at RHIC is highlighted as the first 
of three lines of evidence for QGP production at RHIC. Together with two 
other convergent lines of evidence, jet quenching and the critical D + Au null 
control, I conclude that the QGP has not only been discovered but that a few 
of its remarkable properties have already been established experimentally. 

2. The Empirical QGP 

The discovery of the gedanken QGP phase of matter in the laboratory re- 
quires an empirical definition of the minimal number of necessary and suffi- 
cient conditions in terms of experimentally accessible observables. My empir- 
ical definition is summarized by the following symbolic equation 

QGP = Pqcd + pQCD + dA . (1) 

Why are three independent lines of evidence needed? The first term, PQCD^ 
stands for a class of observables that provide information about its bulk thermo- 
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dynamic equation of state. The equation of state characterizes its long wave- 
length nonperturbative thermodynamic properties briefly reviewed in section 1 . 

The second term, pQCD, stands for class of observables that provide direct 
evidence about its short wavelength dynamics predicted by perturbative QCD. 
The QCD plasma differs quahtative from familiar abeUan QED plasmas due 
to its unique non-Abelian color field dynamics. The radiative energy loss of 
energetic short wavelength partons was predicted to lead to striking quenching 
patterns [13]-[16] of moderate and high px hadrons. The high RHIC cm energy 
of 200 AGeV insures that pt ~ 10 — 20 GeV jet production rates are large 
enough to measure via a wide array of inclusive and correlation observables. 
These hard partons serve as effective "external" tomographic probes of the the 
QGP and test its pQCD c^romo-dynamics. Jets play the analogous role of 
neutrinos that probe the physics of stellar cores, while hadrons play the role of 
photons that probe the corona of the fireball. 

Below RHIC energies, the pQCD line of evidence could not be fully de- 
veloped because the jet rates decrease too rapidly with energy. However, even 
more importantly, at the lower < 4 GeV available the effects of initial state 
nuclear dynamics and the final state hadronic dynamics could not be com- 
pletely deconvoluted from the final spectra. This is the key point that I will re- 
peatedly emphasize which differentiates the observables at the SPS and RHIC 
energies. The necessity to test for the same compHcations at RHIC is what 
gives rise to the third term in Eq.(l). 

The third term, denoted by dA, stands for control experiments that can 
clearly differentiate between alternative nuclear dependences specific to initial 
state partonic wavef unctions as well as the production mechanisms. The con- 
trol differential, dA, is critical at any energy because the QCD plasma must 
first be created from pure kinetic energy ! There has been no hot QGP in the 
universe (except in cosmic ray collisions) since the last drop condensed into 
hadrons about 13 billion years ago. Cold crystalline quark matter may lurk in 
the cores of neutron stars, but the transient hot QGP must be "materialized" 
in the lab. The "matter" arises from decoherence of virtual quantum chromo 
fluctuations in the initial wavefuntions of high energy nuclei. 

At ultra-relativistic energies, these virtual fluctuations are frozen out due 
to time dilation into what has been called a Color Glass Condensate (CGC) 
[17]. The CGC is the high density generalization of the Bjorken-Feynman 
dilute parton model. At high field strengths, the non-linear interactions of 
virtual quantum color fluctuations are predicted to limit the very small Bjorken 
xbj ^ Fourier components. The saturation property of CGC is related to 
unitarity constraints and determines the maximal entropy that can be produced 
in A A at a given ^/s as also pointed out by EKRT [18]. 

The dA control is needed to characterize to what extent these nonlinear 
initial state physics effects can be differentiated from effects due to final state 
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interactions in the QGP matter that forms from it. At RHIC, the best exper- 
imental handle on the dA term happens to be the study of -D + ^ reactions. 
In such light-heavy ion reactions, the initial state CGC physics can be isolated 
because the produced QGP, if any, is too tenuous. 

Why don't I add more terms in Eq.(l)? In fact, each term stands for many 
independent components, as I elaborate below. For the three required terms 
in Eq.(l) the published experimental evidence is now overwhelming and con- 
clusive. Four independent experiments have converged to complementary very 
high quality data sets. 

The three terms in Eq.(l) are necessary and sufficient for establishing that 
a discovery has been made of a uniquely different form of strongly interacting 
QGP. After discussing the three lines of evidence, I will elaborate on why 
I believe that direct photons, J/i), HBT, or other interesting observables do 
not need to be added to Eq.(l). Those observable provide valuable additional 
constraints on the combined and convoluted properties of the initial state, the 
QGP, AND the dense hadronic matter into which it condenses. However, the 
deconvolution of the initial and hadronic final effects has already proven to be 
very difficult at SPS energies and will continue to be at RHIC. 

To avoid misunderstanding, the discovery of the QGP does not mean that its 
physical properties are now understood. In fact, it only signals the beginning 
of a long and well focused direction of research. The history of the neutron star 
discovery offers an instructive analogy. In 1934 Baade and Zwicky proposed 
the theoretical existence of neutron stars soon after Chadwick discovered free 
neutrons. Thirty years later in 1967 Hewish and Bell observed the first few pul- 
sars when suitable radio telescopes could finally be constructed. An amusing 
anecdote is that they actually agonized for a time about whether LGM (little 
green men) were sending them encrypted messages from the cosmos. T. Gold 
in 1968 (as D + Au did at RHIC in 2003) put the debate to rest. Gold pro- 
posed that radiative energy loss of a magnetized neutron star would cause a 
predictable spin down. Later precision measurements confirmed this. Seventy 
years after its proposal, neutron star research still remains a very active exper- 
imental and theoretical direction of physics. Current interest has focused on 
possible color field super-conductivity [19] recently predicted in the very high 
jiB sector of the QCD phase diagram, beyond the boundaries [6] of Fig.2b. 

The critical D + Au control experiments in 2003 could have found that the 
Au + Au QGP observables were strongly distorted by the possible initial CGC 
state that created it. This would have certainly foiled Eq.(I). The search for 
the bulk QGP phase of matter would then have had to await higher energies 
and densities at LHC or for a better understanding of how to deconvolute that 
initial state physics. The large positive signatures in similar p + Pb control 
experiments at SPS showed in fact initial effects strongly distort key observ- 
ables. At SPS the physics of high px Cronin enhancement and p + A ^ J/ijj 
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suppression remain the important open problems. In contrast, at RHIC ener- 
gies, the absence of jet quenching at midrapidity and the "return of the jeti" 
correlations in dA = D + Au provided the check-mate completion of Eq.(l). 

As emphasized by McLerran[17] in these proceedings, the D + Au control 
at RHIC at high rapidities does in fact produce a positive signature for new 
initial state physics. In those kinematic ranges xbj < 0.001, dA fails as a null 
control for QGP, but may signal the onset of nonlinear CGC initial conditions. 
In this lecture, I concentrate on the midrapidity region, xbj > 0.01, where 
Eq.(l) was conclusively satisfied. 

3. Pqcd and Bulk Collective Flow 

The identification of a new form of "bulk matter" requires the observation 
of novel and uniquely different collective properties from ones seen before. 
This requirement is the first term in Eq.(l). In heavy ion collisions, the pri- 
mary observables of bulk collectivity are the radial, azimuthal and longitudinal 
flow patterns of hundreds or now thousands of produced hadrons. Stocker, 
Greiner, and collaborators were the first to predict [9-11, 20, 21] distinctive 
"side splash and squeeze-out" collective flow patterns in nuclear collisions. 
The different types of collective flows are conveniently quantified in terms of 
the first few azimuthal Fourier components [22], Vn{y, PT,Np, h), of centrality 
selected triple differential inclusive distribution of hadrons, h. The centrality 
or impact parameter range is usually specified by a range of associated mul- 
tiphcities, from which the average number of participating nucleons, Np, can 
be deduced. The azimuthal angle of the hadrons are measured relative to a 
globally determined estimate for the colhsion reaction plane angle ^{M). The 
"directed" vi and "elliptic" V2 flow components [12, 22, 23]-[3I] are readily 
identified from azimuthal dependence 

dNhjNp) dNhjNp) 1 ^ „ . , . , 

d^ = ^^2^^' + 2.,(y,p^,iV„^)cos<^ 

+ 2v2{y,PT,Np,h)cos2(j) + ---) . (2) 

The "radial flow" component, "1" , is identified [32] from the hadron mass 
dependence of the blue shifted transverse momentum spectra 



dNh{Np 



2 - exp[-mhCosh{pj_- f3{y))/Tf] , (3) 
't 



dydp, 

where m?i(sinh(/9j_), cosh(pj_)) = (p^, \Jmj^ + p\_) and is the mean 
collective transverse flow rapidity at y. 

Figure (3) shows the striking bulk collectivity elliptic flow signature of QGP 
formation at RHIC. Unlike at SPS and lower energies, the observed large ellip- 
tic deformation ((H- 2v2)/ (1 — '2,V2) ~ 1.5) of the final transverse momentum 
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Figure 3. First line of evidence: Bulk collective flow is the barometric signature of QGP 
production. Left figure combines STAR [26]-[29] and PHENIX [30, 54] measurements of the 
azimuthal elliptic flow {v2{pt)) of ■K,K,p,A. in Au+Au at 200 AGeV. The predicted hydro- 
dynamic flow pattern from [33]-[37] agrees well with observations in the bulk pT < 1 GeV 
domain. Right figure from [24] shows V2 scaled to the initial elliptic spatial anisotropy, e, as a 
function of the charge particle density per unit transverse area. The bulk hydrodynamic limit is 
only attained at RHIC. 



distribution agrees for the first time with non-viscous hydrodynamic predic- 
tions [33]-[43] at least up to about p^- ~ f GeV/c. However, the right panel 
shows that when the local rapidity density per unit area [23, 24] drops below 
the values achieved at RHIC ~ 30/fm^, then the elliptic flow (scaled by the 
initial spatial ellipticity, e = ((y^ — x'^)/{y'^ + x^))) falls below the perfect 
fluid hydrodynamic predictions. We will discuss in more detail the origin of 
the large discrepancy at SPS energies in the next section. 

The most impressive feature in Fig. (3) is the agreement of the observed 
hadron mass dependence of the elliptic flow pattern for all hadron species, 
7r,i^, p, A, with the hydrodynamic predictions below 1 GeV/c. This is the 
QGP fingerprint that shows that there is a common bulk collective azimuthally 
asymmetric flow velocity field, u^{t, r, </>). 

The flow velocity and temperature fields of a perfect (non-viscous) fluid 
obeys the hydrodynamic equations: 

d,{[eQCD{T{x)) + Pqcd{T{x))]u^^{x)u-'{x) - g'^'' Pqcd{T{x))} = , 

(4) 

where T{x) is the local temperature field, Pqcd{T) is the QGP equation of 
state, and eqcoiT) = {TdP/dT — P)qcd is the local proper energy density. 
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The above equations apply in the rapidity window \y\ < 1, where the baryon 
chemical potential can be neglected. Eq.(4) provides the barometric connection 
between the observed flow velocity and the sought after Pqcd- 
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Figure 4. Left figure combines STAR and NA49 data [26] and shows that the directed side- 
wards flow, vi{y), is correlated over 8 units of rapidity at RHIC. At SPS collectivity is domi- 
nated by the overlapping fragmentation regions while at RHIC the nearly identical directed flow 
of in the fragmentation regions is shifted to \y\ > 2. Right figure shows the pseudo rapidity 
dependence of elliptic from PHOBOS [31]. 



The long range nature of collective flow has also been conclusively estab- 
lished by STAR and PHOBOS seen in Fig.(4). The sidewards flow is anti- 
correlated over 8 units of rapidity! In addition, its azimuthal orientation was 
shown to coincide with the azimuthal direction of the largest axis of elliptic de- 
formation at y = 0. This provides an important test of the overall consistency 
of the hydrodynamic origin of flow. The rapidity dependence of the elliptic 
flow in Fig. (4) also shows the long range nature of bulk collectivity. 

Why is V2 more emphasized than vi or radial flow as a signature of QGP 
formation? The primary reason is that elliptic flow is generated mainly during 
the highest density phase of the evolution before the initial geometric spatial 
asymmetry of the plasma disappears. It comes from the azimuthal dependence 
of the pressure gradients, which can be studied by varying the centrality of the 
events [22]. Detailed parton transport [44] and hydrodynamic [39] calculations 
show that most of the V2 at RHIC is produced before 3 fm/c and that ellip- 
tic flow is relatively insensitive to the late stage dissipative expansion of the 
hadronic phase. In contrast, radial flow has been observed at all energies [32] 
and has been shown to be very sensitive to late time "pion wind" radial pressure 
gradients [45], which continue to blow after the QGP condenses into hadronic 
resonances. 
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The observation of near ideal V2 fluid collectivity as predicted with the 
Pqcd together with vi{y) and other consistency checks (c.c.) conclusively 
estabhsh the first term in Eq.(l) : 

Pqcd = V2{pt;t^,K,p,A) +vi{y) + c.c. . (5) 

Preliminary Quark Matter 2004 analysis of H, flow are consistent with the 
predicted V2{pT,y = 0,M,h) and this information is lumped into the c.c. 
terms in Eq.(5). Other data which provide consistency checks of the hydrody- 
namic explanation of collective flow include the observed tt, K, p radial flow 
data [32] for px < 2 GeV. In addition, predicted statistical thermodynamic dis- 
tributions [46] of final hadron yields agree remarkably well with RHIC data. 
Had hadro-chemistry failed at RHIC, then a large question mark would have 
remained about bulk equiUbration in the QGP phase. 

4. QGP Precursors at SPS and Dissipative Collectivity 

It is important to point out, that no detailed 3+ ID hydrodynamic calcula- 
tion [41]-[43] has yet been able to reproduce the rapid decrease of V2{\r)\ > 1) 
observed by PHOBOS in Fig. (4). This discrepancy is due, in my opinion, to 
the onset of hadronic dissipation effects as the comoving density decreases 
with increasing y. From the right panel of Fig.(3), we see that as a decrease 
of the local transverse density from midrapidity RHIC conditions leads to an 
increasing deviation from the perfect fluid limit. The initial density was also 
observed to decrease at RHIC as |?/| increases [47]. Therefore, from SPS data 
alone, we should have expect deviations from the perfect fluid limit away from 
the midrapidity region. It would be interesting to superpose the PHOBOS data 
on top of the NA49 systematics. 

To elaborate on this point, Fig.5 shows CERES data [48] on V2{pt) at SPS 
energy ^/s = 17 AGeV. In agreement with the right panel of Fig. (3), the 
CERES data falls well below hydrodynamic predictions. At even lower en- 
ergies, AGS and BEVALAC, the V2 even becomes negative and this "squeeze 
out" of plane [II] is now well understood in terms of non-equilibrium BUU 
nuclear transport theory [25, 49]. In order to account for the smallness of V2 at 
SPS, hydrodynamics has to be frozen out at unphysically high Tf k,Tc = 160 
MeV. However, the observed radial flow rules out this simple fix. 

The discrepancy of V2 and hydro at SPS energies can be traced to the im- 
portant contribution of the dissipative final hadronic state interactions. The 
hadronic fluid is far from ideal. In approaches [45, 39] that combine perfect 
fluid QGP hydrodynamics with non-equilibrium hadronic transport dynamics, 
the importance of dissipative hadron dynamics at SPS was clearly demon- 
strated. The problem is that the QGP at lower initial densities condenses on 
a faster time scale cannot take advantage of the of the spatial asymmetry to 
generate large V2. The subsequent dissipative hadronic fluid is very inefficient 
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Figure 5. Evidence for dissipative collective flow below RHIC energies. Left: Non equi- 
librium BUU nuclear transport theory [49, 25] can explain the observed elliptic squeeze-out 
(negative V2) collectivity below 4 AGeV. Right: CERES [48] data on elliptic flow at SPS is 
well below hydrodynamic predictions with freeze-out T/ — 120 MeV required to reproduce 
the single inclusive radial flow. Early freeze-out with T/ — 160 MeV, simulating effects of 
dissipation, is needed to reproduce the data. 



in exploiting spatial asymmetry. A factor of two reduction of the initial QGP 
density, therefore, leads to a significant systematic bias of the V2 barometer, not 
only at SPS but also at high |y| at RHIC. Current hadronic transport theory is 
not yet accurate enough to re-calibrate the barometer away from mid-rapdities 
at RHIC. 

In light of the above discussion, the smallness of dissipative corrections in 
the central regions of RHIC is even more surprising. At mid-rapidities, the lack 
of substantial dissipation in the QGP phase is in itself remarkable. Calculations 
based on parton transport theory [44] predict large deviations from the ideal 
non-viscous hydrodynamic limit. Instead, the data show that the QGP is almost 
a perfect fluid. A Navier Stokes analysis [40] is consistent with [44] and shows 
that the viscosity of the QGP must be about ten times less than expected if 
the QGP were a weakly interactive Debye screened plasma. This unexpected 
feature of the QGP must be due to nonperturbative and hence strong coupling 
physics that persists to at least 3Tc. 

4.1 The Minimal Viscosity of the QGP 

One intriguing theoretical possibility being explored in the literature [50] 
is that the shear viscosity, r], in the strongly coupled QGP may saturate at 
a universal super-string bound, r]/a = l/47r. This conjectured duality be- 
tween string theory and QCD may help to explain also the 20% deviation 
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of Pqcd{T) from the ideal Stefan Boltzmann limit. The discovery of nearly 
perfect fluid flow of long wavelength modes with pr < 1 GeV at RHIC is 
certain to fuel more interest in this direction. 

1 propose that a simpler physical explanation of the lower bound on viscos- 
ity follows from the uncertainty principle, as derived in Eq.(3.3) of Ref. [51]. 
Standard kinetic theory derivation of shear viscosity leads to r] = {p{p)X)/3 
where p is the proper density, (p) is the average total momentum, and A is 
the momentum degradation transport mean free path. The uncertainty prin- 
ciple implies that quanta with average momentum components (p) caimot be 
localized to better than Ax ~ Therefore the momentum degradation 

mean free path cannot be defined more accurately than A > l/{p). For an 
ultra-relativistic system, the entropy density is cr 4p, therefore 

which is within 5% of the string theory bound. It is the consequence of the 
universality of the Heisenberg uncertainty principle. Surprisingly, the QGP 
found at RHIC saturates this uncertainty bound and the data clearly rule out 
the order of magnitude larger predictions based on pQCD [51]. See again [44]. 
The long wavelength modes in the QGP are as maximally coupled as ^ = 1 
allows. 

It is "shear" good luck that the mid rapidity initial conditions at RHIC are 
dense enough to essentially eliminate the dilution of elliptic flow due to the 
imperfect hadron fluid formed after the spatial asymmetry vanishes. (Recall 
that {r]/a)H ^ {Tc/Ty/""^ > 1 for T < [51]). At lower energies or 
higher rapidities this good luck runs out and the mixture of near perfect QGP 
and imperfect hadronic fluid dynamics reduces the elliptic flow. 



5. pQCD and Jet Quenching 

In addition to the breakdown of perfect fluid collectivity at high rapidity 
seen in Fig.(4), Fig.(3) clearly shows that hydrodynamics also breaks down at 
very short wavelengths and high transverse momenta, px > 2 GeV. Instead 
of continuing to rise with pT, the elliptic asymmetry stops growing and the 
difference between baryon vs meson V2 reverses sign! Between 2 < pT < 5 
GeV the baryon ^^(pt) exceeds the meson V2^{pt) by approximately 3/2. For 
such short wavelength components of the QGP, local equilibrium simply can- 
not be maintained due the fundamental asymptotic freedom property of QCD. 
I return to the baryon dominated transition region 1 < pr < 5 GeV in a later 
section since this involves interesting but as yet uncertain non-equilibrium non- 
perturbative processes. In this section I concentrate on the px > 2 GeV meson 
observables that can be readily understood in terms QGP modified pQCD dy- 
namics [13, 14]. 
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The quantitative study of short wavelength partonic pQCD dynamics fo- 
cuses on the rare high pT power law tails that extend far beyond the typical 
(long wavelength) scales p < 3T ~ 1 GeV of the bulk QGP. The second major 
discovery at RHIC is that the non-equilibrium power law high px jet distribu- 
tions remain power law like but are strongly quenched [5 2] -[60]. Furthermore, 
the quenching pattern has a distinct centrality, pT, azimuthal angle, and hadron 
flavor dependence that can be used to test the underlying dynamics in many 
independent ways. Below RHIC energies, the initial state Cronin enhancement 



10 



i' 0-10% central; 
I Pb+Pti@\t= 17.3 GeV 
Btaisi~j4=3''»GeV 
V Au+Au@\pi,|=130GeV 
( Au+Au@^=200 GeV 




0.1 



0.01 



10 

PrlGeV/c) 



dN^/dy=200-350 
WA98 7[° (1 7.4 AGeV) 
dN^/dy=800-1 200 
PHENIX 71° (130 AGeV) 
dN^/dy=2000-3500 
▼ PHENIX n" (200 AGeV) 
* STAR (200 AGeV) 




10 
Pt [GeV] 



100 



Figure 6. Jet Quenching at RHIC. Left [61] shows the jet quenching pattern of n discovered 
by PHENIX [52, 53] at RHIC compared to previous observation of high pr enhancement at 
ISR and SPS energies. The nuclear modification factor Raa — dNAA /TAA{b^do'pp measures 
the deviation of AA spectra from factorized pQCD. Right shows predictions [64] of the ^/s and 
PT dependence from SPS, RHIC, LHC based on the GLV theory [65] of radiative energy loss. 



of moderately high px tails was observed in central Pb + Pb reactions at the 
SPS. At the ISR a reduced Cronin enhancement in a + a reactions was seen. 
In contrast, at RHIC a large suppression, by a factor of 4-5, was discovered in 
central Au + Au that extends beyond 10 GeV for vr'^. 

Jet quenching in yl + A was proposed in [15, 16] as a way to study the dense 
matter produced at RHIC energies. As noted before, the pQCD jet production 
rates finally become large enough to measure yields up to high px > 10 GeV. 
Order of magnitude suppression effects were predicted based on simple esti- 
mates of induced gluon radiative energy loss. Ordinary, elastic energy loss [62] 
was known by that time to be too small to lead to significant attenuation. 

As reviewed in [13, 14] refinements in the theory since then have opened 
the possibility of using the observed jet quenching pattern as a tomographic 
tool [63] to probe the parton densities in a QGP. The right panel shows a recent 
jet tomographic analysis [64] of the PHENIX tt^ data [52, 53] based on the 
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GLV opacity formalism [65]. Vitev and I concluded from Fig.6b that the ini- 
tial gluon rapidity density required to account for the observed jet quenching 
pattern must be dNg/dy ~ 1000 ± 200. 

This jet tomographic measure of the initial dNg/dy is in remarkable agree- 
ment with three other independent sources: (1) the initial entropy deduced via 
the Bjorken formula from the measured multiplicity, (2) the initial condition of 
the QGP required in hydrodynamics to produce the observed elliptic flow, and 
(3) the estimate of the maximum gluon rapidity density bound from the CGC 
gluon saturated initial condition [18]. 

These four independent measures makes it possible to estimate the maximal 
initial energy density in central collisions 

eo = e(r ~ 1/po) « 4^ ^ ^ 20^ ^ 100 x sa (7) 
vrit^ dy fm 

where po ~ Qsat ~ 1.0 — 1.4 GeV is the mean transverse momentum of the 
initial produced gluons from the incident saturated virtual nuclear CGC fields 
[17, 18]. This scale controls the formation time fi/po ~ 0.2 fm/c of the initially 
out-of-equilibrium (mostly gluonic) QGP. The success of the hydrodynamics 
requires that local equilibrium be achieved on a fast proper time scale Tgg ~ 
(1 - 3) /po < 0.6 fm/c. The temperature at that time is T{Teq) ~ (eo/(l - 3) x 
12)1/'' ^ 2Tc. 

In HUING model[66], the mini-jet cutoff is po = 2 - 2.2 GeV limits the 
number of mini-jets well below 1000. The inferred opacity of the QGP is much 
higher and consistent with the CGC and EKRT estimates. 

5.1 Raa{pt) and Single Jet Tomography 

In order to illustrate the ideas behind jet tomography, I will simphfy the dis- 
cussion here to a schematic form. See [13] for details. The fractional radiative 
energy loss of a high energy parton in an expanding QGP is proportional to the 
position weighed fine integral over color charge density p{x±, r) 

AEgLv/E^C2k{E) dTTp{x^{T),T) , (8) 

Jo 

where C2 is the color Casimir of the jet parton and k{E) is a slowly varying 
function of the jet energy [13]. The azimuthal angle sensitive escape time, 
L{(j)), depends on the initial production point and direction of propagation rel- 
ative to the elliptic flow axis of the QGP [67]. For a longitudinal expanding 
QGP, isentropic perfect fluid flow imphes that tp(t) « {l/A±)dNg/dy is 
fixed by the initial gluon rapidity density. In this case 
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Therefore, gluon jets are expected to lose about 9/4 more energy than quarks. 
In addition, since the produced gluon density scales as the number of wounded 
participating nucleons at a given impact parameter, (9) predicts a particular 
centrahty and azimuthal dependence of the energy loss. Detailed numerical 
studies show that the actual GLV energy loss can account qualitatively for the 
saturation of V2{pt > 2) [67], the unexpected pr independence [64] of the 
quenching pattern, the centrality dependence [54] of the suppression factor [68, 
69], and the rapidity dependence [74] of Raa{v^Pt) [42, 69] for pions and 
high > 6 GeV inclusive charged hadrons at RHIC. 

To further illustrate qualitatively how (9) influences the quench pattern, 
consider a simplified initial jet distribution rate, (fN/(fpo = cp^^ where 
n ~ 7. In applications these are of course calculated numerically. After 
passing through the QGP, the final jet pr = Pf = Po(l ~ e)- The average 
over fluctuations constrains (e) = AE{(f))/E. The quenched jet distribution is 
(fN/(fpT = {(fpo/(fpf)(fN/(fpQ « (1 - e((^))"-2 cpj". The calculated 
hadron inclusive distribution is obtained by folding the quench jet distribution 
over the fragmentation function D{z = PhjPj^Q'^ = P/)- However, in this 
illustrative example the fragmentation function dependence drops out, and the 
nuclear modification factor, Raa^Pk-, 0; ^part) = d^hi^) / dNh{0) reduces to 

The average over e takes into account fluctuations of the radiative energy loss. 
The resulting Raa is independent of pT in this approximation. This was also 
found in the detailed numerical work in [64]. In central collisions, Npart ~ 2 A 
and L{(l)) ^ (L), and the magnitude of quenching is fixed by ec oc dNg/dy. 
The centrality and azimuthal dependence for non-central collisions follows 
without additional calculations. 

5.2 Iaa and Di-Jet Tomography 

Measurements of near side and away side azimuthal angle correlations of 
di-jet fragments provide the opportunity to probe the evolution of the QGP 
color charge density in even more detail. Fig.(7) show the discovery [57- 
59] of mono-jet production [15] in central collisions at RHIC. In peripheral 
collisions, the distribution dN/ dAcp of the azimuthal distribution of ~ 2 
GeV hadrons relative to a tagged ~ 4 GeV leading jet fragment shows the 
same near side and away side back-to-back jetty correlations as measured in 
p + p. This is a direct proof that the kinematic range studied tests the physics 
of pQCD binary parton collision processes. For central collisions, on the other 
hand, away side jet correlations are almost completely suppressed. 
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Figure 7. Monojets at RHIC from STAR [58, 59, 57]. Strongly correlated back-to-back di-jet 
production in pp and peripheral AuAu left side is compared to mono-jet production discovered 
in central AuAu. 

The quantitative measure of the nuclear modification of di-jet correlations in 
A + B reactions at a given ^/s is given by a formidable multi- variable function 

CABiyi,PTl,(l>l,y2,PT2,4>2;b,^b,hl,h2) , (11) 

where (yi , Pti , ) is the trigger particle of flavor hi and the {y2,PT2,4'2) 
is an associated particle of flavor /12 for collisions at an impact parameter b 
with a collective flow axis, the reaction plane, fixing the azimuthal angle $5. 
Obviously Cab is a very powerful microscope to study the modification of 
short wavelength correlations in the strongly interacting QGP. 

The published data are as yet limited to yi ~ 2/2 ~ 0, broad px cuts: pTi > 
4 GeV and p-p ~ 2 GeV, two bins of (pi — (p2, and of course averaged over 
The measured modification of di-jet correlations is obtained by subtracting out 
the correlations due to bulk elliptic flow via the di-jet measure 

Iaa = d{cl)i-(l)2){N{<t>i-(l)2) 

-Nb{1 + 2v2{pi)v2{p2)cos{2{(I)i-cI)2))} , (12) 

where the number of triggered pairs N{4>i — 4>2) is normalized relative to the 
expected number based on p + p measurements in the same [A_, A+] rela- 
tive azimuthal angle range. Wang [70, 71] has analyzed the centrality, Npart, 
dependence of Iaa as well as Raa and showed that both can be understood 
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from the same pQCD energy loss formalism. This provides another critical 
consistency check of jet tomography at RHIC. 

Additional preliminary data from STAR presented by K. Filimonov at Quark 
Matter 2004 showed the first direct evidence that back-to-back dijet quenching 
has a distinctive dependence on the azimuthal orientation of the jets relative to 
the reaction plane as expected from the obvious generaUzation of Eq.(lO) 

Iaa{A(I> = TT, ^b) ~ ({(1 - hLi{ro,PT)){l - hL2{ro, -PT))r-') (13) 

where rb is the initial transverse production point of the approximately back-to- 
back dijet moving in directions pr relative to the reaction plane ^i,. Here ki, the 
effective fractional energy loss per unit length in the QGP produced at impact 
parameter b. The high pxi trigger naturally biases fb to be near the surface 
and pt is biased toward the outward normal direction. This means that on the 
average, Li <C L2, and the away side fragments should be strongly suppressed 
in the most central collisions, while the near side fragment correlations should 
be similar to that seen in pp. However, in non-central minimum biased events, 
Eq.(13) naturally predicts the away side fragments are less quenched when the 
trigger hadron lies in the reaction plane than perpendicular to it, as observed 
by STAR. 

5.3 The Empirical pQCD Line of Evidence 

Single and dijet data for pions above 3 GeV and protons above 6 GeV pro- 
vide conclusive evidence that the QGP matter is partially opaque to short wave- 
length probes with a quenching pattern as predicted by the pQCD radiative 
energy loss: 

pQCD = RauAu{PT, 4>, Npart) + IAuAu{<Pi - (t>2] b, $5) . (14) 

With the vastly increased statistics from the current RHIC RUN 4, the tests of 
consistency of the theory will be further extended to px ~ 20 GeV, and Iaa 
and Caa will attain ever greater resolving power. In addition, heavy quark 
tomography [76-78] will provide new tests of the theory. 

6. The dAu Control 

Only one year ago [79] the interpretation of high px suppression was under 
intense debate because it was not yet clear how much of the quenching was 
due to initial state saturation (shadowing) of the gluon distributions and how 
much due to jet quenching discussed in the previous section. There was only 
one way to find out - eliminate the QGP final state interactions by substituting 
a Deuterium beam for one of the two heavy nuclei. In fact, it was long ago 
anticipated [16] that such a control test would be needed to isolate the unknown 
nuclear gluon shadowing contribution to the A-i-A quench pattern. In addition 
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Figure 8. The dA control: PHENIX [72] 7r° and STAR [73] data compare RdAu to 
RauAu- These and BRAHMS [74] and PHOBOS [75] data prove that jet quenching in Au+Au 
must be due to final state interactions. Curves for tt'^ show predictions from [64] for AuAu and 
from [82] DAu. The curves for DAu show the interplay between different gluon shadow 
parameterizations (EKS, none, HIJING) and Cronin enhancement and are similar to predictions 
in [80-82]. In lower panel, the unquenching of charged hadrons is also seen in D + Au relative 
to Au + Au at high pr- 



D+Au was required to test predictions of possible initial state Cronin multiple 
interactions [80-82, 84]. In contrast, one model of CGC [83] predicted a 30% 
suppression in central D+Au. 

The data [72-75] conclusively rule out large initial shadowing as the cause 
of the xbj > 0.01 quenching in Au+Au and establish the empirical control 
analog of Eq.(14) 



dA = Rdau{pt, 4>, Npart) + IdAu{4>i - 4>2,b) 



(15) 
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6.1 The "Return of the Jeti" 

The IdAu measurement from STAR [73] is the check mate! The return of 




Figure 9. The dA "Return of the Jeti": Dijet fragment azimuthal correlations from STAR [73] 
in DAu are unquenched relative to the mono jet correlation observed in central AuAu. 

baclc-to-baclc jet correlation in D + Au to the level observed in pp is seen in 
Fig.9. The data rule out CGC gluon fusion models that predict mono-jets [86] 
correlations in the xbj > 0.01 region. These D + Au data support the conclu- 
sion [70, 71] that the observed jet quenching in AuAu is due to parton energy 
loss. 

Another independent check of the strikingly different nature of the nuclear 
matter created in Au + Au versus D + Au is provided by the width of the 
away side correlation function [87]. The transport properties of cold nuclear 
matter extracted from the Cronin enhancement effect [82] have only a small 
effect on the measured dijet acoplanarity. Preliminary PHENIX data presented 
by J. Rak at Quark Matter 2004 confirm the qualitative similarity ofp + p and 
D+Au for pt > 2 GeV, but at the same time demonstrate a strong quantifiable 
increase in acoplanarity in central Au + Au consistent with multiple semi-hard 
scattering [88] in a dense QGP 

7. Conclusions 

The three lines of evidence have converged from the four RHIC experi- 
ments. The empirical QGP that has been found at RHIC via the combination 
of Eqs.(l,5,14,15). This QGP is, however, not the weakly interacting [7], color- 
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dielectric "wQGP", that we have searched for. Because of its near perfect fluid 
long wavelength properties, it must be very strongly coupled at least up to sev- 
eral times Tc- Symbolically, we should denote ^ the empirical QGP at RHIC 
by "sQGP" to emphasize its special properties. The sQGP is not only a near 
perfect fluid but it also retains part of its QCD asymptotic freedom character 
through its highly suppressed, but power law, short wavelength spectrum. 

In summary, the sQGP found RHIC was seen through the following three 
convergent lines of evidence 

sQGP = PqcD+pQCD + dA 

= {v2{pt]t^, K,p, K) +vi{y) + c.c.} 

+{RAuAu{PT, (I>, Npart) + lAuAu{,(f>l - 4>2;b, 
+{RdAu{PT, ^, Npart) + lDAu{(t>l - <t>2,h)} ■ (16) 

Other surprising properties are already known. The anomalous pt = 2— 5 GeV 
baryon/meson ratios [27, 89] , noted in connection with Fig. 3 and also seen in- 
directly in Fig. 8, already point to unexpected novel baryon number physics. 
Current speculations center around possible gluonic baryon junction dynam- 
ics [92], and possible multi-quark quark coalescence mechanisms [93]. The 
baryon number transport properties in the sQGP will certainly teach us new 
physics. 

The experimental task of mapping out the novel properties of sQGP has only 
begun. It is important to concentrate, however, on those observables which are 
least distorted or "polluted" by uninteresting hadronic final interactions. For 
example, the severity of the HBT puzzle depends on which hadronic transport 
model is used. Pure hydrodynamics with late freeze-out times fails badly to re- 
produce final state soft pion correlations. Hybrid hydro+RQMD hadronic cas- 
cade does somewhat better [90], but there is at least one transport model [91] 
that reproduced the data. Non-equiUbrium hadron resonance transport dynam- 
ics are unfortunately still not well enough understood at any energy to allow 
definitive conclusions to be drawn. 

Thermal direct photons have yet to be measured, but it is akeady known that 
the pre-equiUbrium [80] and hadronic final state contributions in the few GeV 
Pt range will produce large backgrounds on top of the thermal component. 
These must be deconvoluted if thermal photons are to serve as a sQGP ther- 
mometer. However, even the theoretical thermal photon rates are still under 
debate [94]. 

The J /ifj suppression discovery at SPS was originally attributed to Debye 
screening of cc in a wQGP paradigm. Recent lattice QCD results now indicate 

'l thank T.D.Lee for discussions and suggesting tlie sQGP designation of the matter discovered at RHIC to 
distinguish it from weakly interacting, Debye screened plasmas, wQGP, which may only exist at tempera- 
tures T > Tc. 
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that heavy quark correlations persist perhaps up to 2Tc. This is another indica- 
tion for the strongly coupled nature of sQCD. The suppression of J/V' in AA 
is also strongly influenced by initial state and final state hadronic (comover) 
effects. It will be no easier to deconvolute these competing effects at RHIC. 
These observables of course need to be measured at RHIC, but one should not 
expect an easy interpretation. 

I believe that the most promising direction of future experiments at RHIC 
will be precision measurements of short wavelength (px > 2 GeV) corre- 
lators Cab illustrated in Eq.(II). These are very powerful six dimensional 
microscopes with four discrete {hi, /i2, A, B) and two continuously adjustable 
geometric (6, experimental knobs in addition to the beam energy ^/s = 
20 — 200 AGeV. One of the important correlators will be that of direct photon 
tagged jets [95]. Another will be open charm and possibly bottom quark to- 
mographic probes. The available experimental knobs have hardly been varied 
yet. Much remains to be done to map out to clarify the properties of the sQGP 
found at RHIC. 
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